CU-TP-952 



How well do domain wall fermions realize chiral symmetry? 
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In the domain wall fermion formulation, chiral symmetry breaking in full QCD is expected to fall exponentially 
with the length of the extra dimension. We measure the chiral symmetry breaking due to a finite extra dimension 
in two ways, which can be affected differently by finite volume and explicit fermion mass. For quenched QCD 
the two methods generally agree, except for the largest extent of the extra dimension, which makes the limit 
uncertain. We have less data for full QCD, but see exponential suppression for the method where we have data. 



1. Introduction 

Lattice QCD with massless domain wall 
fermions (including fermion loop effects) [p]]0 is 
expected to have the S\J L {N f ) ® SU R (Nf) chiral 
symmetry of the continuum when the extent of 
the extra dimension, L s , becomes infinite. For 
simulations, where the volume is finite and parti- 
cles are not strictly massless, reliable techniques 
are needed to quantify the symmetry breaking for 
finite L s . Such techniques are needed to see the 
expected exp(— aL s ) dependence of chiral break- 
ing for full QCD and determine if this is also the 
case for the quenched theory. 

Here we report results from two techniques for 
measuring chiral symmetry breaking due to finite 
L s ; the first uses the pion mass and the second 
the axial Ward identity. At zero temperature, the 
pion mass is governed by the axial Ward identity. 
However, in simulations with finite volume and 
with finite quark masses, it is important to check 
the agreement between these approaches. 

The axial Ward identity is the origin of the 
Gell-Mann-Oakes-Renner (GMOR) relation, dis- 
cussed previously for domain wall fermions in || . 
The fermion action of M is used, with the mod- 
ifications of Q. Some details on the numerical 
methods are in || . See || for a general review on 
domain wall fermions and references. 
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2. m ros for domain wall fermions 

If the dominant effect of finite L s is to produce 
an extra contribution, m ros , to the total quark 
mass, then one would expect 

ml = c (V) + c x {V) *(m f + m rcs ) + • • • (1) 

where V is the space-time volume and it is ex- 
pected that cq(V) — > as V — > oo. For finite 

(m 2 ) 

volume, a result we call m^cs , can be found 
from m\{mf — 0)/ci(V), which is m res when 
co (TO - 0. 

Using the flavor non-singlet axial current in (^] 
we integrate the axial Ward-Takahashi identity 
to get 

(<7o<?o) = w/Xtt + A J 5 (2) 

where the pseudoscalar susceptibility is (no sum 
on a) 



Xt, 



A" 



A" 



(3) 



the additional contribution from chiral mixing 
due to finite i, is 



AJ 5 = 2j2{35(x-Ls/2) m 
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(4) 



and q x are four-dimensional fermion fields defined 
from the appropriate right- and left-handed fields 
at the boundaries of the extra dimension. 

For large volumes in the chirally broken phase, 
the pseudoscalar susceptibility is expected to be- 
have as 



Xtv = <x_i/(m/ + m rcs ) + a + 0(m f + : 



(5) 
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The first term again says that, for large volumes, 
the pion is massless at m/ = — m rcs , while ao 
gives the contribution due to the massive modes. 
Clearly, the pion pole contribution only domi- 
nates for large enough volumes and small enough 
to/ + m ros . 

jS;(x,L s /2), a pseudoscalar density located 
midway between the domain walls, also has a pole 
contribution, whose coefficient is suppressed by 
propagation from L s /2 to the boundaries. Since 
Xtt and AJ5 both have a pole at to/ = — m rcs and 
when the pole terms dominate (Q) {qoQo) is finite, 
we can write in general 

AJ 5 = TOresX, + 60 + 0(m f + TO res ). (6) 

We define m les ^ GMOR ' by simultaneously fit- 
ting to the form 

(9o9o) = (mf + m ICS )x-K + b a . (7) 

and Xtt as given in (||). For a given L s , this is 
a four parameter fit for a_i, ao, &o and mres** ^ • 
Note that only if bo/x-n is small, can we get a reli- 
able estimate for mif s MOR ^ from (qoqo) /x-K~ m f- 
For full QCD, both mi^ s M ° and bo should ap- 
proach zero exponentially in L Sl since both in- 
volve propagation from L s /2 to the walls. 

3. m res in quenched QCD 

We first find mrcs"'' for the the — 5.7, too = 
1.65, 8 3 x 32 quenched domain wall spectrum 
study we reported last year 0. For quenched 
QCD, the observed zero mode contribution to 
(<Zo<7o) is small for to/ > 0.02, so we restrict our 
attention this mass range here. Figure pi shows 

TOres for our (3 — 5.7, to = 1.65, 8 3 x 32 simula- 
tions from a correlated, linear fit of rri^ to valence 
masses 0.02, 0.06, 0.10, with errors from jacknif- 
ing. The L s — 32 and 48 values are the same 
within errors, making the large L s limit seem non- 
zero. For L s ~ 24 the result for a 16 3 x 32 lattice 
is also shown, revealing that finite volume effects 
are noticeable. 

Figure [l] also shows m^ B MORS> from a corre- 
lated fit to (0) and f|. The 16 3 GMOR point 
is on top of the 8 3 point in the plot. The fits 
TO rcs = 0.059(14) exp[-0.052(10)L s ] and b = 
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Figure 1. Comparison of m Tes ( m ^ and 
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Figure 2. (qq) jx-n and valence m rc / GMOR ' on 
8 3 x 4,/3 = 5.2, to = 1.9, to/ = 0.02 lattices. 

-0.0035(11) exp[-0.051(13)L s ] (not shown) are 
consistent with AJ5 vanishing in the L s — > 00 
limit. 

The agreement between the two methods for 
L s < 32 is reasonable and only occurs since 60 
is included in the fits. m^ s MOR ^ is volume inde- 
pendent for L s = 24, while TOrS is not. The 
discrepancy for L s = 48 may be due to finite vol- 
ume, but needs further study. 

4. mres in TV/ = 2 QCD 

For full QCD, we have done extensive simu- 
lations with the Wilson gauge action and do- 
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Table 1 

Valence m ros comparison on 8 3 x 32 lattices 





m f 


(mi) (GMOR.) , 

mros mics -b 


Wilson gauge action, Q — 5.325, mo = 1.9 


24 


0.02 
0.06 


0.0622(9) 0.058(2) 0.0047(3) 
0.0645(6) 0.059(2) 0.0046(2) 




Iwasaki 


gauge action, 8 = 1.9, mo = 1-9 


24 


0.02 


0.0401(5) 0.038(2) 0.0028(2) 




Iwasaki 


gauge action, = 2.0, to = 1.9 


24 
48 


0.02 
0.06 
0.02 


0.0158(9) 0.015(2) 0.0013(3) 
0.019(1) 0.017(3) 0.0015(4) 
0.0073(16) 0.011(2) 0.0010(2) 



main wall fermions on 8 3 x 4 volumes for sev- 
eral values of L s with 8 — 5.2, mo = 1.9 and 
to/ = 0.02. For these lattices, which are in the 
low temperature phase, we show (qq) jx-w in Fig- 
ure 0. An exponential fit, yielding (qq) /xtt ~ 
0.02+0.082(3) exp[-0.027(2)L s ] for 16 < L s < 40 
with x 2 /Ndof = 2.76/2, is also shown. 

We have also done uncorrelated fits, using 
valence masses between 0.02 and 0.14, to ex- 
tract m[es M ° R ^ and b , since for the dynami- 
cal simulations there is not enough data to re- 
solve the covariance matrix. (For the quenched 
case there was little difference between the cor- 
related and uncorrelated fits.) All fits have 



(GMOR) 



IS 



A do f = 4 and x /AW < 1. m tl 
also shown in Figure @ and the dashed line fit is 



???.,. 



(GMOR) 



0.17(2) exp[-0.026(6)L s ] for 10 < 
L s < 40 with x / N dof = 0.35/4. b is not 
shown, but also fits the exponential form bo = 
-0.0100(16) exp[-0.0147(67)L s ] with x 2 /A dof = 
0.20/4 over the same range in L s . 

In Table [J we compare the two methods of ex- 
tracting TO r es using valence spectrum data from 
Nf = 2, 8 3 x 32 scale setting calculations ||. All 



data was fit for 0.02 < m 



(val) 



< 0.1. For the 



GMOR fit, A dof = 2 and x 2 /N do f < 1 for all 
fits. We note that the two methods agree within 
statistics. 

We can also calculate m rcs both ways but with 
data as a function of the dynamical mass. Since 
there are only two dynamical masses, both meth- 
ods are unconstrained so the errors quoted come 
from naive extrapolation. The results are sum- 
marized in Table @. 



Table 2 

Dynamical m rcs comparison on 8 3 x 32 lattices 
3 m (ml) m (GMOR) -6n 

fj fibres ''<<rcs uq 



m = 1.9, L a = 24 



5.325 0.059(2) 0.053(7) 0.004(1) 
2.0 0.013(2) 0.014(5) 0.0011(9) 

5. Conclusions 



We have gotten good agreement between 
and m^ s MO for a wide range of quenched and 
dynamical simulations by including the non-pole 
terms in the susceptibilities. From our current 
data, 77irS M ° fl ' appears less volume dependent. 
For the quenched simulations at L s = 48, the two 
methods do not agree, possibly as a result of finite 
volume effects. This case warrants further study. 

Whether chiral symmetry is fully restored for 
quenched simulations in the L s — > oo limit of do- 
main wall fermions is still an open question. For 
Nt 



2 QCD, m[g MOR) 



falls exponentially, even 



at quite strong coupling. The rate of chiral sym- 
metry restoration is very slow leaving much room 
for improvement. 

All numerical calculations were performed on 
the 400 Gflops QCDSP machine at Columbia and 
the 6 Gflops QCDSP machine at Ohio State. 
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